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Abstract. Surface waves on a stationary flow of water are considered, in a linear 
model that includes the surface tension of the fluid. The resulting gravity-capillary 
waves experience a rich array of horizon effects when propagating against the flow. 
In some cases three horizons (points where the group velocity of the wave reverses) 
exist for waves with a single laboratory frequency. Some of these effects are familiar 
in fluid mechanics under the name of wave blocking, but other aspects, in particular 
waves with negative co-moving frequency and the Hawking effect, were overlooked 
until surface waves were investigated as examples of analogue gravity [Schiitzhold R 
and Unruh W G 2002 Phys. Rev. D 66 044019]. A comprehensive presentation of the 
various horizon effects for gravity-capillary waves is given, with emphasis on the deep 
water/short wavelength case kh^ 1 where many analytical results can be derived. A 
similarity of the state space of the waves to that of a thermodynamic system is pointed 
out. 



PACS numbers: 04.70.Dy, 04.70.-s, 05.45.-a, 05.70.Fh, 47.35.-i, 92.05.Bc, 92.10.Hm 
1. Introduction 

The interest in black-hole analogues has been mainly driven by the intriguing possibility 
of observing Hawking radiation in the laboratory [U El [3l HI [5l [6] . In addition to the 
experimental challenges, this pursuit has important theory implications because of the 
well-known weakness in the derivation of the Hawking effect for real black holes [71 18] . 
Hawking's semi-classical calculation [9] is based on a consideration of fields that are 
assumed to have no appreciable gravitational effect compared to the black hole, but 
the derivation contradicts this assumption because the fields attain arbitrarily high 
frequencies (and therefore energies) at the horizon. This so-called trans-Planckian 
problem reveals the lack of any proper understanding of quantum-gravitational effects. 
Unfortunately, the question of whether black holes really radiate does not seem to be 
amenable to experimental investigation because such radiation would be completely 
swamped by the cosmic microwave background. In black-hole analogues the trans- 
Planckian problem is avoided by dispersion that limits the blue-shifting of waves at 
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horizons [in]. The Hawking effect has been found to persist even in the presence of 
dispersion [TOl [TTl [121 HSl Ell [HI HSl |T7], although in general numerical simulations 
must be resorted to in order to verify this. The lack of a good understanding of 
the Hawking effect with arbitrary dispersion is an important and presumably solvable 
problem, whereas a proper understanding of the Hawking effect in real black holes 
requires unknown physics to deal with the infinite blue-shifting. Analogue systems 
therefore provide the opportunity to understand the Hawking effect using known physics 
and to verify it experimentally; neither of these tasks can be achieved with real black 
holes. 

The introduction of dispersion into the Hawking effect opens up a variety of 
new possibilities. Indeed, one of the lessons of analogue systems is that the physics 
of horizons, including the Hawking effect, has a breadth and richness that is not 
immediately apparent from the case of real black holes. One can consider white-hole 
and black-hole horizons [HI [6], that both give rise to the Hawking effect; horizons allow 
the blue-shifting and red-shifting of probe waves [6]; depending on the dispersion, the 
probe waves can bounce back at horizons, or go straight through, or bounce back and 
forth a number of times (see below); two horizons can communicate leading to runaway 
quantum Hawking radiation or a classical instability [HI [HI [5l [20] . 




Figure 1. Experimental white-hole horizon in hydrodynamics. The surface waves 
propagate against the flow up to the point where the flow speed matches the group 
velocity of the waves. The flow speed is higher on left than on the right because of the 
slope on the bottom of the tank. 

Perhaps the least exotic black-hole analogue that has been proposed is that of 
surface waves on a moving fluid. Schiitzhold and Unruh [IB] showed that long- 
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wavelength surface waves in a shallow moving fluid obey the Klein-Gordon equation 
in a curved space-time geometry. By varying the flow speed so that it exceeds the 
wave group velocity in some region, one produces a horizon. Close to the horizon 
the long-wavelength assumption breaks down and the detailed behaviour is governed 
by the dispersion; the wave evolution is no longer described by an effective space- 
time metric and there is therefore no trans-Planckian problem of infinite blue-shifting. 
Further results for surface waves from an analogue-gravity perspective were given in [21] . 
Recently, we performed experiments to investigate the behaviour of surface waves at a 
white-hole horizon, using water waves in a wave tank with a counter-flow [22]. The 
speed of the counter-flow was varied along the tank by the presence of a sloping region 
on the bottom. Figure [T] shows how waves propagating against the flow are blocked 
at the point where the counter-current reaches the group velocity of the wave; this 
blocking line is the white-hole horizon. Video of the incoming wave was used to produce 
the space-time diagram in Figure |2] (note that the directions of the flow and the wave 
are reversed in the plot of Figure |2] compared to Figure [T]). The linear features with 
positive slope in the space-time diagram are the evolution of the wave crests and troughs 
in time — the world-lines of the crests and troughs. The inverse dx/dt of the slope of 
the world-lines gives the speed of the crests and troughs, the phase velocity u/k of the 
wave since the phase is J {kdx — wdt); the lines curve upward as the wave reaches the 
white-hole horizon, showing a decrease in the phase velocity. The wavelength in the 
laboratory is revealed by drawing a horizontal line in the diagram and measuring the 
distance between two crests; because of the curving upward of the crest world-lines, the 
wavelength is seen to decrease as the horizon is approached. This is the characteristic 
wavelength-shortening ( "blue- shifting" ) of waves at a white- hole horizon, also observed 
in an optical analogue [6J. The inset in Figure |2] shows the behaviour of rays at a 
white-hole horizon where there is no dispersion; the rays stick at the horizon which 
corresponds to an infinite blue-shifting of the waves (the trans-Planckian problem). It 
is the dispersion of the surface waves that limits the amount of blue-shifting at the 
horizon; in this paper we will discuss some of the rich set of possible behaviours. 

Not surprisingly, the blocking of waves by a counter-flow that exceeds the wave 
group velocity is a phenomenon that is well-known in the fluid-mechanics community [231 
[35| |25| l29l [30] . The blue-shifting of waves at the blocking line (white-hole horizon) is 
also well-known and has been investigated experimentally [271 [2H1 129] . There is however 
another possible process in the interaction of a wave with a counter-flow, one that does 
not seem to have been considered by the fluid-mechanics community, even though it 
is a possibility clearly visible in the dispersion relation [THl 122]: this process is the 
Hawking effect. In a stationary flow, where the flow speed at each point in the tank is 
constant in time but varies from point to point, the frequency in the laboratory frame 
of a surface wave on the flow is conserved. As we have seen, the wavelength in the 
laboratory is not conserved, so a white-hole horizon converts the wavelength of a wave 
propagating against the flow to a different wavelength while conserving its frequency in 
the laboratory. Figure [2] is an example of such a wave; it is right-moving against the 




Figure 2. Space-time diagram of a plane wave encountering a white-hole horizon. In 
this diagram the wave propagates to the right on a flow moving to the left. Grey values 
denote the water height and show clearly the world-lines of the crests and troughs. The 
green dotted line shows the initial slope of the world-line of the incoming wave, which 
is inversely proportional to the phase velocity. The red line is the world-line of a crest, 
the slope of which increases as the wave reaches the horizon. The phase velocity thus 
decreases at the horizon, leading to a decrease of the wavelength, a blue-shifting that is 
limited by dispersion. The inset shows the behaviour of rays at a non-dispersive white- 
hole horizon, where there would be infinite blue-shifting of waves (rays originating on 
both sides of the horizon are shown). 



flow (ingoing against the flow) with positive (angular) frequency u and positive wave 
number k and the horizon increases k while keeping u fixed. An important quantity is 
the frequency of the wave in a frame co-moving with the flow; this co-moving frequency 
is not conserved, but the ingoing wave and the blue-shifted wave with higher k both 
have positive co-moving frequencies [TH1122]- It turns out that there is often a solution 
of the dispersion relation for the flxed positive input frequency u that has a negative k 
and a negative co-moving frequency [181 122] . The laboratory frequency u of the input 
wave must be conserved in the wave evolution but when there is a wave with negative 
co-moving frequency at same value of u, there exists the possibility that it could be 
generated in the interaction of the input wave with the counter-flow. These waves 
with negative co-moving frequency would be produced in addition to the blue-shifted 
waves with positive co- moving frequency; this process is the Hawking effect, which is at 
root a classical effect, although with extraordinary quantum implications |9l|38]. The 
strangest feature of the Hawking effect is that it is an ampliflcation of the ingoing wave, 
an extraction of energy from the flow (or from whatever background creates the horizon). 
The famous quantum Hawking radiation is a rather straightforward consequence of the 
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classical Hawking effect when the fields underlying the waves are quantized (the role of 
input waves is then played by the quantum vacuum and the energy extraction in the 
Hawking process allows the spontaneous creation of field quanta) [9] . 

In our experiments described above [22] we observed indications of waves with 
negative co-moving frequency. Numerical simulations indicated that the Hawking effect 
would be unobservable in the regimes covered in the experiments [22], but it is not 
clear how well the wave evolution was described by the theoretical model used in the 
simulations. It appears that our paper was the first experimental search for the Hawking 
effect with surface waves, and further experiments are needed. We hope through this 
paper to increase awareness of the Hawking effect among researchers in fluid mechanics 
by showing how it has been overlooked in the existing literature on wave blocking. 
Complementary to this, we hope to show those familiar with horizon physics some of 
the rich horizon effects that occur in models of wave blocking that are used in fluid 
mechanics. Such surprising connections between apparently disparate areas of physics 
(and engineering) are often a source of inspiration to both sides. 

It is important to stress that the model of surface waves considered in this paper 
is a linear one: the dynamics is given by a linear wave equation, but with a nonlinear 
dispersion relation. As is well known, surface waves behave non-linearly if the wave 
amplitude is large enough, and this presents a potential problem since the shortening of 
the wavelength at a horizon necessarily leads to an increase in the wave amplitude. The 
limiting of blue-shifting by dispersion also limits the amplitude increase at blocking lines 
and, for the linear model to be applicable, it is important that the wave amplitude not 
get too large. The maximum amplitude attained during the interaction with a counter- 
flow can to some extent be controlled through the amplitude of the incident wave, and 
dissipation (see next paragraph) also helps to limit the amplitude. The interaction of 
nonlinear waves with a counter-flow is a highly complex process that is not discussed 
here (see for example [31]). 

Also neglected in the model used below is the viscosity of the fluid, and this 
limits the accuracy of the results when wavelengths are blue-shifted down to very small 
(capillary) wavelengths since viscosity quickly damps such waves. But it is easy to 
make allowance for this latter drawback of the model by bearing in mind the limited 
propagation length of capillary waves. 

2. Surface waves on a stationary flow 

The subjects of black holes and water waves have a historical connection through 
the figure of Pierre-Simon de Laplace. In his "Exposition du Systeme du Monde" 
in 1796, Laplace famously introduced the term etoile sombre (dark star) to denote an 
object whose gravitational field is strong enough to prevent light from escaping (the 
same concept had been described in 1783 by the Reverend John Michell in a letter to 
Cavendish) |32j. Laplace is also well known for having derived (in 1775) a dispersion 
relation for surface waves on water |33j. When the water has a background flow (which 
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does not vary with depth) and the effect of surface tension is included, the dispersion 
relation takes the form widely used in fluid mechanics to describe waves on moving 
water [231 [HI ESI El] : 

(u; - Ukf = (^gk + tanh(fc/i). (1) 

Here u is the (angular) frequency of the wave in the laboratory frame and k the wave- 
number; U is the speed of the flow, h the depth of the fluid, g the gravitational 
constant, p the fluid density and 7 the surface tension. For water, p = 1000 kg m~^ 
and 7 = 0.073 Nm-^ Equation g is a one-dimensional dispersion relation suitable 
for water-tank experiments. When the flow is stationary {U independent of time, 
but gradually varying in space), w is a constant but k varies with spatial position x. 
Waves described by ([T]) are a consequence of gravity (g) and surface tension (7) and are 
called gravity-capillary waves. For small k the gravity term dominates, which we call 
the gravity regime, whereas for large k surface tension dominates, giving the capillary 
regime. The pure-gravity case corresponds to 7 = 0. 

The quantity uj — Uk is the frequency in a frame co-moving with the fluid. Hence 
the positive, respectively negative, square roots of ([T]) 




uj-Uk = ±^l {gk + ^k^jtanh{kh) (2) 

correspond to positive, respectively negative, co-moving frequencies. As described in 
the Introduction, the Hawking effect is the generation of a wave on the negative branch 
of Q from a wave on the positive branch, through interaction with a counter-flow. It 
is a remarkable fact that in the extensive fluid-mechanics literature on the waves ([T]), 
including the pure-gravity case 7 = 0, there seems never to have been any consideration 
of the possibility of conversion of waves from the positive to the negative branch of ^ 
through the blocking effect [23l|3a[35l[25l[26l|36l|37l[27l[28ll^ If such a 

conversion process had been investigated, the Hawking effect would presumably have 
been (re-) discovered in fluid mechanics. 

Solutions of the dispersion relation ([T| are usually represented graphically. We 
obtain from ^ 

uj = Uk± ^ (^gk+^k^^ tanh(fc/i) . (3) 

Figure |3] plots both branches of the right-hand side of ^ as functions of k, for a 
fixed value of [/ < 0; figure (a) shows the pure-gravity case (7 = 0) and figure (b) 
shows the full dispersion relation with surface tension included. The positive branch 
of ([3]), corresponding to positive co-moving frequency, is shown in green while the 
negative branch, corresponding to negative co-moving frequency, is shown in blue. The 
intersection of these curves with a given horizontal line, such as the red line in the figures, 
gives the possible waves for the frequency u given by that line. In a stationary flow 
CO is conserved but the plots of the right-hand side of ^ change with spatial position 
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as U{x) changes and one can trace the evolution of a given solution by following its 
intersection point with the horizontal line of fixed u. As one traces the evolution of the 
intersection point, the changing group velocity dw/dfc is the slope of the tangent to the 
curve at the point of intersection with the horizontal line. 





Figure 3. Dispersion relation for surface waves propagating on a water flow with a 
given velocity U < 0; (a) shows the pure gravity case (7 = 0) and (b) shows the gravity- 
capillary case. The green curves are the positive branch of the right-hand side of ([s]) 
(positive co-moving frequency) and the blue curves are the negative branch (negative 
co-moving frequency). Intersections of these curves with the horizontal red line (value 
of Lu) show possible waves at that lo. If U becomes more negative the green and blue 
curves rotate clockwise about the origin. Local maxima and minima of the green and 
blue curves show the possibility of blocking waves by means of appropriate velocity 
profiles U{x) (see Figure |4]). A water depth h = 0.4m is assumed. 
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In the pure-gravity case (Figure [3] (a)) there are at most four real roots of the 
dispersion relation (fc/, fc^, kji and fc/^). The solution fc/ is a right-moving wave in the 
laboratory, having positive phase and group velocities, propagating against a left-moving 
counter-flow U <Q. If the flow speed increases as the wave moves to the right (f/ 

becomes more negative) the green and blue curves tip over clockwise about the origin 
(as in Figure |4]), so the root ki increases — the wave is blue-shifted. When the wave 
reaches a point where the flow speed has increased to make the roots ki and ks coalesce 
at a local maximum of the green curve, the group velocity of the wave is zero — it has 
been blocked at a white-hole horizon. The wave has been stopped by a negative "group 
acceleration" that is still non-zero at the blocking point so the group velocity decreases 
to negative values; the wave moves back to the left in the laboratory on the fcs root 
of the dispersion relation, back into the region where the counter-flow is slower than 
the blocking speed and where the wave previously had wave number ki. The ingoing 
wave ki has thus been blue-shifted to ks by the white-hole horizon; the blue-shifted 
wave ks has positive phase velocity, so its crests move to the right in the laboratory, 
but it has negative group velocity. The third real root of the dispersion relation with 
positive co-moving frequency (green curve) in Figure [s] (a) is k^] this is simply a wave 
propagating in the same direction as the flow, to the left with negative phase and group 
velocities. The solution k^i is rather trivial and is of no interest for horizon effects. 

The root ku in Figure [s] (a) has negative co-moving frequency (blue curve) and is 
of great interest for horizon effects; the conversion of some of the input wave ki to ku 
is the Hawking effect. Since it has a negative wave number, the wave kn has a negative 
phase velocity in the laboratory; its crests move backwards relative to the direction of 
the ingoing wave, in contrast to ki and fc^, which both have positive phase velocity. It 
is essential to understand that the three waves fc/, ks and kn are all propagating to 
the right relative to the fluid, even though in the laboratory k^ has a phase velocity 
pointing left and both ks and ku have a group velocity pointing left {k,R is the root 
corresponding to a wave moving to the left relative to the fluid). Unlike the blue- 
shifting of ki to ksi the existence of conversion from kj to ku cannot be deduced from 
dispersion plots, which only reveal it as a possibility. The amount of conversion of 
kj to kn depends on the details of the dispersion and the velocity profile U{x). In 
simple cases involving limited dispersion the Hawking effect is determined by the slope 
d[/(x)/dx of U at the horizon (this slope is the analogue of the surface gravity of a 
black hole, the acceleration due to gravity at the horizon), but for general dispersion 
and velocity profiles no analytical formula for the size of the effect has been found and 
one must resort to numerical simulations of the wave evolution. One aspect of the 
challenge to find a good intuitive understanding of the Hawking effect is apparent from 
the description of the horizon given above: this was taken as the point where the flow 
speed matched the group velocity of the blue-shifting wave. But the phase velocity of 
this wave is greater than its group velocity, so one can have a group-velocity horizon but 
no phase- velocity horizon. On the other hand, one can have both a group- and a phase- 
velocity horizon, with in principle an arbitrary distance between these two horizons and 
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a completely different value of (iU{x)/(ix at each horizon. The size of the Hawking 
effect is influenced by these and other factorsj|] There is also the further possibility of 
having the maximum flow speed close to but less than that required for a group- velocity 
horizon. In this case one would expect some wave tunneling into the blue-shifted root 
/cb, and perhaps also into (tunneling of surface waves has been studied in [iU]). 
Numerical simulations indicate that this method of generating kn without a group- 
velocity horizon is mathematically possible for a steep enough velocity profile, but it 
should not be possible in practice [22]. In the experiments reported in [22] waves with 
negative phase velocity were observed even in the absence of a white-hole group- velocity 
horizon, but, as stated in the Introduction, the origin of those waves is not clear. 

The conversion of ki to ks discussed above is well known in fluid mechanics, under 
the name of wave blocking [231IME5l[23[2SlESlEZll2Zl[2Hl[2SEniE3. ^ superposition 
of the ki and ks waves has been shown to be describable by an Airy interference 
pattern [HI |26l |36l [37]. In contrast, the root kn in Figure |3]^a) has been largely 
neglected by the fluid-mechanics community. Although the graphical representation of 
the dispersion relation is standard in fluid mechanics, very few authors [251 Il2l HSl HI] 
plot the negative-fc part, in either the pure-gravity or gravity-capillary cases, and the 
conversion of ki to kn appears not to have been considered. 

Turning to the full gravity-capillary case. Figure [3](b) shows (for flxed [/ < 0) how 
the surface tension 7 changes the dispersion relation at large wave numbers compared 
to the pure-gravity case in Figure [s] (a). For the value of U plotted, the positive co- 
moving frequency curve (green) has a local minimum as well as a local maximum, and 
this is also the case for the negative co- moving frequency curve (blue), although the 
local minimum of the latter curve always occurs at negative laboratory frequency uj. 
If U becomes more negative the curves tip over clockwise about the origin, so that for 
larger counter-flow speed there exist roots with negative co-moving frequencies at the u 
shown by the red line. Each local maximum or minimum of the green and blue curves 
reveals the possibility of reversing the group velocity of a wave with an appropriate 
velocity proflle, as in the discussion of the local maximum in the pure-gravity case 
(Figure [3] (a)). The local maximum of the green curve in Figure [s] (b) allows the wave 
blocking and blue shifting of an incident right-moving wave as in the pure-gravity case. 
We refer to this blocking line as the white horizon (from white hole). But after the 
white horizon has reversed the group velocity of the incident wave so that it now moves 
to the left, this blue-shifting wave encounters another blocking line because of the local 
minimum of the green curve. We refer to this second blocking line for the blue-shifting 
wave as the blue horizon. At the blue horizon the group velocity reverses once more to 
become positive so the wave moves to the right again towards the white horizon. This 
time the (still blue-shifting) wave goes right through the white horizon, so overall the 
incident wave undergoes a double bounce. Figure |4] shows the graphical solution of the 
dispersion relation at four values \Ui\ < \U2\ < jf/sl < \Ui\ of a velocity proflle U{x) < 

I An extension of the standard analytical results to the dispersive case was given in PT; but only for 
two specific velocity profiles, both of which gave a group- and a phase-velocity horizon. 
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Figure 4. Graphical solution and numerical ray solution for gravity-capillary waves 
on a water counter-flow. The green/blue curves refer to positive/negative co-moving 
frequency. Four values of the velocity profile U{x) < are shown in the dispersion 
plots (upper) and the a;-positions where the profile takes these four values are shown 
in the ray solutions (lower). The values are (in ms^^) Ui = —0.17, U2 = —0.2035, 
Us — —0.2536 and U4 = —0.275. We use a hyperbolic tangent function tanhax 
to describe the variation of the velocity profile, with a typical length a = 0.5 m. 
The conserved frequency uj of the waves (red horizontal line in the dispersion plots) 
corresponds to a period T = 0.7 s. A water depth h = 0.4 m is assumed. 
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and also a numerical solution of the ray equations for the incident right-moving wave 
(green curve in the lower figure). The x-positions where the velocity profile takes the 
four values used in the dispersion plots are shown in the ray plot. Rays move at the 
group velocity and so the wave blocking is clear from the ray plots. Also shown is the 
ray solution for the wave with negative co-moving frequency (blue curve); this wave is 
initially left-moving in the laboratory but its group velocity is reversed at a blocking 
line we refer to as the negative horizon (this horizon does not exist in the pure-gravity 
case — see [22] )• Comparing the dispersion and ray plots in Figure |4] (and ignoring the 
co-propagating wave discussed above) one can see how the increasing counter-flow 
speed as x increases gives, successively, one k root in the profile where \U\ < |t/2|, three 
roots in the region where U lies between U2 and f/3, five roots in the region where U lies 
between f/3 and [4, and three roots in the region where \U\ > \U4\. One can also see 
how these roots relate to the ray behaviour. (See |22] for ray plots in the pure-gravity 
case.) 

Figure [5] shows an example where the frequency u and counter- flow profile U{x) < 
are such that only a wave with positive co-moving frequency exists. The wave again 
displays the double-bouncing behaviour. Figure |6] shows a numerical solution for a wave 
packet centered on the ray in Figure [5j This simulation was obtained by solving the 
scalar wave equation describing the surface wave on the counter-flow; it was shown 
by Schiitzhold and Unruh that this equation takes the form of the Klein-Gordon 
equation in a curved space-time with added higher-order dispersion, in this case the 
dispersion of gravity-capillary waves: 



The method of numerically solving equations of the form (|4]), for essentially arbitrary 
dispersion, is described in [10]; further examples of this kind of numerical solution for 
surface waves appear in [TH] and [22] • In the wave-packet simulation in Figure [6| the 
continuous blue-shifting that accompanies the double bounce is apparent. Because of 
the spread of frequencies in the wave-packet there is some leakage of the initial wave 
through the white horizon (first bounce) as well a spreading and separation of frequency 
components at the blue horizon (second bounce). For extensive numerical simulations 
of gravity-capillary waves in the presence of a current, see HI]. 

In practice the blue-shifting of incident gravity waves into the capillary regime, as 
in Figure |6j will be limited by viscosity, which is not included in the model we have been 
discussing. As a consequence, the highly blue-shifted waves produced at the blue horizon 
will dissipate rapidly. An experimental investigation of these effects was first performed 
by Badulin et al. Pl]. These authors observed the initial blocking of waves in the gravity 
regime (white horizon) and the subsequent conversion into waves in the capillary regime 
(blue horizon) which then propagated through the original blocking line and vanished 
through viscous damping [34j . Gravity-capillary waves on a counter-flow have also been 
studied experimentally by Klinke and Long [36], who produced space-time diagrams of 
the wave evolution. Trulsen and Mei [3S1 133] give a theoretical treatment that includes 
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Figure 5. Graphical solution and numerical ray solution for gravity-capillary wave 
with positive co-moving frequency. The conserved frequency uj corresponds to a period 
T = 0.47s. For this frequency there is no wave with negative co-moving frequency in 
the range \Ui\ < \U2\ < {Usl of counter-flow speeds used. The three values of the flow 
velocity shown are (in ms^^) Ui = -0.1804, U2 = -0.1876, U3 = -0.1991. A water 
depth h = 0.4 m is assumed. 

numerical simulations. A recent theoretical survey is given by Huang 

We have so far discussed only certain features of the dispersion relation ([T]). In 
the next Section we classify in more detail how the presence or absence of the various 
horizons, and their positions, depend on the conserved frequency u and velocity profile 
U{x). 

3. Results for wavelengths less than the water depth 

The influence of the water depth h in the dispersion relation ([T]) disappears when 
\k\h ^ 1, so that tanhkh ~ ±1 (depending on the sign of k). This is the case of 




Figure 6. Wave packet simulation. The packet is centred on the ray 



in Figure [Sj 
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wavelengths short compared to the water depth and it gives the polynomial dispersion 
relation 



that is easier to handle analytically than the orig inal The waves considered in 
Figure [3] to Figure [6] are in fact very well described by the deep water/short wavelength 
dispersion relation ([s]), as were the waves studied in the experiments [22] and [M] . 

At the blocking lines or horizons discussed in the last Section, the group velocity 
vanishes and the dispersion curve u{k) has a local extremum. Three possible horizons 
for gravity-capillary waves were identified: the white, blue and negative horizons in 
Figure [3|b). At each horizon two real roots of the dispersion relation coalesce into one 
double root and then disappear: in the terminology of dynamical systems it is a saddle- 
node or tangent bifurcation [H]. The order parameter of the bifurcation is the wave 
number whereas the two control parameters are the velocity U (the "external field") 
and the frequency uj (the "internal parameter" ) . Following the approach in [H] , we find 
the horizons by looking for a double root /c2 of the cubic dispersion relation ([s]): 



where ki is the remaining simple root. Comparing coefficients of /c in ([5]) and (|6]) we 
obtain expressions for ^2 and fci, as follows. The comparison of coefficients gives three 
equations for the two unknowns ki and A;2; two of these equations are solved for ki and 
/c2, and the the third equation is then a constraint relating ki to k2- In the case of 
positive wave numbers (plus sign in ([s])) this procedure gives 







(6) 




(positive wave numbers) 



(7) 




with the constraint 

kik\ = - — (positive wave numbers). 

7 

The constraint ([s]), for both sign possibilities in ([7|), leads to 
' ^5 gU^ -fUJ^U^ 157wf/2 Qg-fU 7^^ 2772^2- 







positive wave 
numbers 



(8) 



(9) 



|_ 4a; pg 2p p pu Ap'^g 

For negative wave numbers (minus sign in (|5])) ^2 and ki are 




(10) 



with the constraint 




(negative wave numbers) 




7 
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that leads to 



P9 



2p 



p 



negative wave 
numbers 



(12) 



The significance of the constraints l\9h and (12) is clear from Figure 3[b). With a 



particular choice of U, local extrema in the dispersion plot (corresponding to double 
roots ^2) occur at values of u determined by this choice of U. For a given U, the 
constraint (|9]) or (12) is a quartic in to whose real roots give all the frequencies at which 
a blocking line (horizon) occurs for this U. Alternatively, upon fixing u the constraint 
gives a quintic in U whose real roots are all the fiow velocities that give a horizon at 
this frequency. 




T(s) 



Figure 7. The flow speeds \U\ at which the white, blue and negative horizons occur, 
as a function of the wave period T (black curves). The curve lying on or close to 
the straight red line is the white horizon, the red line showing the pure-gravity white 
horizon. The curve approaching the asymptote \U\ = \Uj\ from below is the blue 
horizon, and the curve approaching this asymptote from above is the negative horizon. 



As in the previous Section, we consider only positive u and negative U, whereas 
k can be positive or negative (this gives no restriction in the horizon effects). The 
white horizon in Figure occurs at a value of k low enough for the influence of the 
surface tension 7 in (§ to be very small. A very good approximation for the white 
horizon, which is exact in the pure-gravity case, is therefore obtained by putting 7 = 0. 
Neglecting 7, the constraint ^ gives 

ioU\U+^)^0. (13) 



Note that the approximation (13) is also obtained in the large U limit, so for large U 
it becomes the exact constraint at the white horizon for gravity-capillary waves. From 
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(13) we recover the relation between U and uj that gives the blocking of gravity waves 
at a white-hole horizon Hi] 



f/» = -f = . (14) 

where T is the period. The corresponding value of the double root k2 is obtained from 
the 7 — )■ limit of ([T]); the lower sign in the first of ([T]) gives the only finite expression: 

^ 9 



Equations (14)-(15) show the exact relationship between the frequency, wave number 
and counter-flow speed at the white horizon for pure-gravity waves. Note from (14) 
that the flow speed at the white horizon is proportional to the conserved period of the 
blocked wave. The straight line \Ug\ versus T is shown in red in Figure [7j Superimposed 
on that red line is a black curve that shows the exact relationship between flow speed 
and period for gravity-capillary waves. In line with the comments above, the red line 
agrees very well with the gravity-capillary case except for small |[/|, and therefore small 
T. The striking feature of the gravity-capillary curve is that it ends at the point labeled 
{\Uc\,Tc)] the white horizon thus does not exist for periods T that are below a critical 
value Tc, or for counter- flows that do not reach a critical speed \Uc\- The existence of 
this threshold can be seen from the dispersion plots in Figure |4j for flow velocity Ui 
there is no local maximum of the green curve, so no frequency u experiences a white 
horizon at this flow velocity; in contrast, the other flow velocities plotted in the figure all 
give a local maximum of the green curve and therefore a white horizon for the frequency 
at this maximum. Similarly, if the period T is too small (frequency uj too large), the 
horizontal red line in Figure |4] will not intersect a local maximum in the dispersion plot 
for any f/, so there can be no white horizon for such periods. 

Let us look in more detail at the critical values (|?7c|, Tc). Figure [s] shows graphically 
the occurrence of the threshold for the white horizon. We see that the disappearance of a 
local maximum in the green curve, as \U\ decreases, is accompanied by the disappearance 
of the local minimum, and at the critical value Uc the two local extrema coalesce to form 
a point of inflection at frequency ujc- This shows that (|f/c|,Tc) is also the threshold for 
the occurrence of the blue horizon, which requires a local minimum in the green curve. 
At (|?7c|,Tc) the two double roots k2 in ([T]) (white and blue horizons) coincide, and in 
fact the same value is taken by the simple root fci, as can also be seen from Figure [8] 
The values (|f/c|, Tc) can be obtained by solving for the point of inflection ^ = = 
in the dispersion relation (|5| (/c > 0); these two equations can be solved for Uc and the 
critical wave number kc-, and Uc then follows from the dispersion realtion. Alternatively, 
the point of inflection is found by demanding that the square-root expression in ([t]) 
vanishes so that all /c2 and ki coincide. The result is 

/ \ 1/4 

T, = 27r(3 + 273)3/4 { —A = 0.425 s (16) 
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Figure 8. Dispersion plots and ray solutions for a wave with period Tc = 2-k jujc- The 
wave with positive co- moving frequency occurs at a point of inflection when U — Uc- 
This means that the corresponding ray (green) has a group velocity that slows to zero 
when \U\ increases to \Uc\', the group velocity does not reverse, however, and the ray 
resumes its propagation into regions of higher \U\. 



and 



\/3 f "yg^ 



Uc = ' , — =-0.178m/s (17) 



with the wave number 

k. = ^r^V'' = 144m- 



(3 + 2^3)1/2 

In using the method ^ of searching for the system parameters at horizons, we noted that 
they correspond to saddle-node or tangent bifurcations in dynamical-systems theory 
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The cusp {\Uc\,Tc) in the \U\ vs T diagram (Figure [7|) is the point where two saddle-node 
hues (horizons) intersect. In the terminology of dynamical systems this corresponds to 
a so-called pitchfork bifurcation |47j . 




Figure 9. The dispersion plot when the counter-flow velocity is such that the local 
minimum of the green curve and the local maximum of the blue curve lie on the k- 
axis {ijj — T — oo). This counter-flow velocity is given by (20) and has the value 
U-^ = -0.231 m/s. 



Since {\Uc\iTc) is also a threshold for the existence of the blue horizon, the curve 
in the \U\ versus T diagram (Figure [?]) relating the period to the flow speed at the blue 
horizon must also end at (|?7c|,Tc); this curve is also shown in Figure [T] and it is seen to 
approach an asymptotic value, labelled ?Xy, as T — oo. This is because T oo means 
a; — and the dispersion plot in Figure [9] shows that the local minimum giving the blue 
horizon occurs at a; = for a finite non-zero U that we call U^. We find the flow velocity 
U^i as follows. The constraint (|9]) relates the values of U and uj at all horizons (double 
roots of the dispersion relation) for waves with positive k. Hence by taking w — )■ in 
(|9| we obtain f/^; this limit of ([9| gives 



0, (19) 
P 



so 



/ \ 1/4 

= ( = -0.231 m/s. (20) 

The wave vector at the local minimum (blue horizon at w = 0) in Figure [9] is found by 
inserting the velocity into the expression for the double root /c2 in ^ and taking the 
upper sign: 

(21) 
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The other local extremum at A; > in Figure [9} the white horizon for a non-zero u given 
by the horizontal red line in the figure, occurs at 

/ \ 1/2 

A;2 = 0.137 (^J , (22) 
and the corresponding u (red line) is 

u = 0.180 (^^^ ) (23) 

where the exact but lengthy numerical coefficients have not been reproduced. The 
simple root ki for U = U-y is zero for cu = 0, while for u given by (23) (red line Figure |9]) 
the simple root is (intersection of red line with green curve at large A; > in Figure |9]) 

/ \ 1/2 

ki = 1.73 (^j . (24) 

These last three results are obtained from the constraint ^ and the expressions Q for 
the double and single roots, with U = U^. 

We see from Figure [9] that in the limit w — )■ (T — )■ oo) the negative horizon also 
occurs at the flow velocity U-y, as well as the blue horizon, and the wave vector at the 
negative horizon is minus that at the blue horizon, —ky. Figure [9] shows that, for waves 
with positive laboratory frequency u, the counter-flow velocity Uy is the threshold for 
the existence of waves with negative co-moving frequency; the threshold flow velocity 
for such waves is the threshold for the negative horizon. It follows that the curve in the 
\U\ versus T diagram (Figure [T]) relating the period to the flow speed at the negative 
horizon must lie above the line \U\ = \U^\ and asymptotically approach this line as 
T ^ oo. This negative-horizon curve is also plotted in Figure [7) it lies above the 
blue-horizon curve but shares with it the asymptote \U\ = \Uy\. Unlike the white- and 
blue-horizon curves in Figure [7| which both end at the cusp {\Uc\,Tc) for small T, the 
negative-horizon curve diverges to |[/| — )■ oo as T — )■ 0. This behaviour of the negative 
horizon is clear from the dispersion plots because as T — )■ (w — t- oo) the flow speed 
\U\ must increase without limit in order for the local maximum in the negative-fc curve 
to reach the horizontal frequency line (see for example Figure |4]). 

The flow velocity Uy that appears as an asymptote in the \U\ versus T diagram 
(Figure [T]) has additional signiflcance in fluid mechanics. Firstly, it is a well-known 
property of gravity-capillary waves on static water {U = 0) that the minimum phase 
velocity of the waves is given (apart from the sign) by the expression (20) for Uy. The 
velocity Uy is also important in the case of shear flows, i.e. velocity proflles that change 
with the fluid depth. It was shown by Caponi et al. [48j that a sufficient condition for 
a shear flow to become spontaneously unstable is for the flow velocity on the surface of 
the fluid to exceed Uy] the instability leads to the generation of gravity-capillary waves 
on the fluid surface |3H] . Another example in shear flows is the appearance of negative- 
energy waves at the interface of two fluid layers, which occurs when the relative velocity 
of the layers exceeds U^; this is related to the famous Kelvin-Helmholtz instability, as 
discussed by Fabrikant and Stepanyants 
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Figure 10. Dispersion plots and ray solutions for waves with period T < T^. The 
Roman numerals I and V refer to counter-flow speeds that lie in the regions labeled 
by these numerals in Figure [7| The period is T = 0.382 s and the flow velocities are 
-0.159 m/s (I) and -0.277 m/s (V). 



Figure [7] allows the classification of the behaviour of gravity-capillary waves on a 
stationary counter-flow. The period T is conserved in the wave evolution so by fixing a 
vertical line in the figure one can distinguish five qualitatively different possibilities for 
waves of a single frequency: 

1. T < Tf.. As one moves into regions in the velocity profile U{x) < with higher 
counter-flow speeds, one moves from region I in Figure [7] into region V. The dispersion 
plots and ray solutions for this case are shown in Figure [lOj Notable is the fact that 
there is no white or blue horizon. 
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2. T = Tc. Here the line of constant T is the vertical green line in Figure [7] that 
separates region I from regions II and III and passes through the cusp point {\Uc\iT^. 
The dispersion and ray plots are shown in Figure [8) as already discussed, this case is 
the threshold for the appearance of the white and blue horizons. 




Figure 11. Dispersion plots and ray solutions for waves with period Tc < T < Ti,. The 
Roman numerals refer to counter-flow speeds that lie in the regions labeled by these 
numerals in Figure]?] The period is T = 0.510 s and the flow velocities are — 0.169 m/s 
(II), -0.194m/s (IV), -0.231 m/s (III) and -0.267m/s (VI). 



3. Tc < T < Tb. The line of constant T lies between the green and brown vertical 
lines in Figure [7] and so passes through region III. Here increasing counter-flow speeds 
takes us from region II to IV to HI to VI in Figure [7j The dispersion plots and ray 
solutions for this case are shown in Figure 11 Here there is a white and blue horizon, 
and the white horizon occurs at a lower counter-flow speed than the negative horizon. 
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Figure 12. Dispersion plots and ray solutions for waves with period T ^ Ti, (see 
Figure [?]). Both the white and negative horizons occur at the same counter-flow 
velocity U — Ub, the point {Ti,, \Ub\) being the intersection of the curves for these 
two horizons in Figure [t] For water, Tt = 0.647 s, Ub — — 0.255 m/s. 



4. T = Th, defined by the brown vertical line in Figure [7) The line T = Tb passes 
through the point where regions III, IV, VII and VI meet, at \U\ = \Ui,\. The significance 
of the point {\Ub\,Tb) is that it is the intersection of the white- horizon curve and the 
negative-horizon curve; this means that for a wave with period Tb the white horizon 
occurs at the same counter-flow speed as the negative horizon. Figure 12 confirms this 
in the dispersion plots and ray solutions. The values {Ub, Tb) must be found numerically 
and for water they are Tb = 0.647 s, Ub = —0.255 m/s. 

5. T > Tb. Here the line of constant T passes through region VII in Figure [7j 
Increasing counter-flow speeds takes us from region II to IV to VII to VI. The dispersion 
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Figure 13. Dispersion plots and ray solutions for waves with period T > Tfc. The 
Roman numerals refer to counter-flow speeds that lie in the regions labeled by these 
numerals in Figure]?] The period is T = 0.692 s and the flow velocities are — 0.192 m/s 
(11), -0.214m/s (IV), -0.260 m/s (VII) and -0.281 m/s (VI). 



plots and ray solutions for this case are shown in Figure 13 Here there is a white and 
blue horizon, and the white horizon occurs at a higher counter-flow speed than the 
negative horizon (compare carefully with case 3 above). 

Trulsen |31j inferred from his results the structure of the (|t/|, T)-diagram of 
Figure [7} but without the negative-horizon curve. He also derived the existence of 
the cusp {\Uc\,Tc) as a triple- root solution of the dispersion relation [S]. 

Recently, we observed (with continuous waves trains) the regions VI, VII and IV 
of Figure [T] as reported in our experimental \U\ versus T diagram [22]. The distinction 
between the regions II and IV was unclear from our data. Our wave-maker was limited 
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to a minimum period of 0.5 s which is higher than Tc = 0.4255 s. Our focus was on 
the conversion into waves with negative co-moving frequency rather than conversion 
into capillary waves in the double-bouncing behaviour described above. Conversion to 
capillary waves is difficult to observe experimentally because of the rapid dissipation of 
capillary waves (we were unaware of the work by Badulin et al. [M])- In addition, we 
used a rather high period (far from Tc) to get long wavelengths of the ingoing waves 
since the waves with negative co-moving frequency should be produced with a drastic 
reduction of the wavelength according to the dispersion relation. We were surprised to 
find indications of waves with negative co-moving frequency even without wave blocking 
(a white horizon). 

Badulin et al. [M] performed experiments with wave packets (three to ten wave 
cycles centred on periods in the range T = 0.33-0.66 s) sent on a counter-fiow with speed 
\U\ between 0.04 and 0.3 m/s over a sloping bottom. Double bouncing of the input waves 
was observed with a strong reduction in both wavelength (from A = 0.2m to 2 mm!) 
and amplitude. These authors presented only one photograph, at T = 0.52 s, of the 
conversion phenomenon and one measurement of the amplitude of waves as a function 
of the position/velocity (T = 0.5 s and dU/dx = 0.1 s^^), but beautiful measurements of 
the effect of the velocity on the wavelength. No results for periods less than Tc = 0.425 s 
were reported. 



4. A thermodynamic analogy 

In the previous Section we summarized the behaviour of gravity-capillary waves (with 
kh ^ 1) on a stationary counter-fiow by means of a diagram in the (|?7|, T)-parameter 
space (Figure [?]). This diagram allows one to visualize the evolution of an incident 
wave of a single frequency (which is conserved), as was illustrated in the five cases in 
Figures 10 8, [TTp2 and 13 In this Section we note a similarity between Figure [7] and 



a phase diagram in thermodynamics, where the horizon lines in Figure [7] are analogous 
to the lines separating different phases (first-order phase transition). In particular, the 
cusp point {\Uc\,Tc) in Figure [T] looks like a critical point (second-order phase transition) 
in a phase diagram. 

Let us explore this thermodynamic analogy a little further. In thermodynamics a 
system is described by an equation of state of the form /(P, V, 9) = 0, where P is the 
pressure, V is the volume and 6 is the temperature. For example, the equation of state 
of a van der Waals gas can be written [47] : 

V'-[nh+—y' + —V- — = ^, (25) 

where n is the number of molecules divided by Avogadro's number, R is the gas constant, 
h relates to the non-zero volume of the molecules and a is a measure of the molecular 
interaction. A familiar property of the van der Waals gas is the existence of a critical 
point in the (P, 0)-phase diagram; this is associated with a fold catastrophe in the 
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surface (25) in (P, V, 9)-space that constitutes the state-space of the gas |17]. Now the 
dispersion relation ([s]), written as 



AttU 



k 



0, 



(26) 



7 7 V" ^ y 7^^ 

describes the state of a wave on a counter-flow as a surface in {U, k, T)-space and this 



surface has the same kind of fold catastrophe as the van der Waals gas — see Figure [14 
The connection between wave blocking and catastrophe theory was inferred a long time 
ago by Peregrine and Smith [36] and more recently by Trulsen [Hj. The white and 
blue horizons appear as curves on the state surface in Figure [14| which converge and 
join at the "critical point" {Uc,kc,Tc) given by (16)-(18). Only the positive-fc part of 



the state surface is shown in Figure [T4J since this contains the point analogous to the 

on to the (f/, T)- 



14 



thermodynamic critical point. By projecting the surface in Figure 
plane one obtains the part of the "phase diagram" Figure [7] containing the white and 
blue horizons, similar to the (P, B)-diagram of the van der Waals gas [17] . 

We can summarize the analogy between the van der Waals gas and waves on 
a counter-flow by the following table, which we stress describes only a qualitative 
relationship: 



Van der Waals Gas 


Wave-Current Interaction 


Volume V 


Wave Number k 


Temperature 


Frequency u 


Pressure P 


Flow Velocity U 


Compressibility |^ 


Susceptibility ^ 


Spinodal Line 


Blocking Line 


Perfect Gas P - 


Pure Advection f/ = f 



The perfect gas is seen to correspond to very large u and U (pure advection U = u/k of 
the surface waves). In addition, the perfect gas is obtained by setting the parameters a 
and b to zero; similarly, pure advection of the surface waves corresponds to setting the 
parameters 7 and g to zero. 



The fold catastrophe in Figure 14 can also be projected to the {U, A;)-plane 
(Figure 15). From the dispersion relation ^ these projections are given by 



Uik) 



k 



9 , 7, 
Y + -k 
k p 



(27) 



with a fixed value of lo (T). Figure 15 shows the curves U{k) for different periods T, for 
both the gravity-capillary case (red lines in (a), (b) and (c)) and the pure-gravity case 



(blue lines in (a), (b) and (c)). Figure 15 ^d) shows U{k) for a range of periods T in the 
gravity-capillary case; these isoperiod curves are analogous to the Andrews isotherms 
for a real gas. 

We have seen from dispersion plots that wave-blocking corresponds to local extrema 
of the function u{k). It follows from the implicit function theorem that the minimum 
(maximum) of U{k) corresponds to the maximum (minimum) of u}{k), and therefore 
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Figure 14. The state of a wave on a counter-flow is a surface in {U, k, T)-space given 



by (26 1. This surface has a fold catastrophe in the positive-fc sector, similar to the 



fold catastrophe of the van der Waals gas [17] . The red curve on the state surface is 
the white horizon and the blue curve is the blue horizon. The horizon curves meet 
at {Uc, kc, Tc), which is analogous to the critical point of the van der Waals gas. The 
white- and blue-horizon curves in Figure [7] are the projection of the curves on the state 
surface to the {U, T)-plane. 



blocking lines are given by the local extrema of U{k). These local extrema dU/dk = 
are the analogues of the spinodal line dP/dV = of the van der Waals gas. Without 
surface tension (blue lines in Figure [l5|^a)-(c)), U{k) has a single minimum that lies on 
the dotted black line in these plots. With surface tension (red lines in Figure [l5|(a)-(c) ) , 
U{k) has in addition a local maximum. At T = Tc (Figure [Tsj^b)) the minimum and the 
maximum of U (k) merge and an inflection point appears; this corresponds to the cusp 
in the ([/, T)-plane (Figure [7|. We can define a mechanical susceptibility Xm = {§if)j'^ 
analogous to the isothermal compressibihty coefficient Xe = —y that diverges at 

the horizons, just as the compressibility of the gas diverges at the spinodal line. 

In the limit of infinite period T — t- oo (w — > 0) we found that the wave is described 
by a fiow speed (20) and wave number (21). We recover this result from the u = case 
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Figure 15. The red curves in Figures (a), (b) and (c) show t/ as a function of k for 
different periods T: Figure (a) is for T > Tc, (b) for T ^ Tc and (c) for T -> oo. 
Figure (d) shows a series of plots of U{k) for different T; the green curve is for T — )■ 
(pure advection). The blue curves in Figures (a)"(c) show the pure-gravity case. The 
black dotted line in all four figures shows, in the pure-gravity case, the locus of the 
local minimum of U{k) as T changes. The horizontal dotted purple line is U^; this is 
the value of U at the local maximum of U{k) for T — )■ oo (Figure (c)). The horizontal 
dotted green line is Uc\ this is the value of U at the point of inflection of U{k) for 
T ^ Tc (Figure (b)). The horizontal dotted orange line is L/g, which is proportional 
to T (see (14)); it gives the value of U at the local minimum in the pure gravity case 
(blue curves in Figures (a)-(c)). 



of (27), which is 

U{k) = - 
This is plotted in Figure 



9 , 7, 
J + -k. 
k p 
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c) and has a maximum ^ 



at wave number 



(2J 



(29) 



where Ic is the capillary length. This wave number in (28) reproduces U^, given by (20). 

The principle of corresponding states implies that the properties of real gas are 
universal functions of the state variables scaled to the critical point. For the Van der 
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Waals gas, it is well known that the equation of state can be written in a universal form 
around the critical point (V^ = 3n6, 6c = gf^R and Pc = 2^)- 

Pr + ^) {3Vr - 1) = 8^,, (30) 



using the scalings {kr = = 3.i'^4 and f/,. = -f-^y^), we find the universal 



where the subscript r means reduced variable (K = ~ W' ~ w)- Similarly, 

using the scalings ( 
dispersion relation: 

{Ur-UrKf = Kil + k^,). (31) 

Whatever the fluid (surface tension, density), its wave-like behavior will be the same 
close to the cusp. The dimensionless form of the constraint ^ becomes: 

UrU^ + hj^ + ulU^ - ^-ulU^ - 6uJrUr - 1 - = 0. (32) 

/ \ 1/4 

One recovers = —a/2 ( ^ j by imposing Ur = 0. By introducing another scaling 
UL = UrOJr = — the constraint reads: 



U'l + \u't + cot (U'l - y f/'J - 6t/; - 1 - '^^cot ] = (33) 
One recovers Ug = —^^ hj imposing Ur = 0. 



5. Conclusions and Perspectives 

We have described the interaction of linear gravity-capillary waves with a counter- 
flow, with emphasis on the various horizon effects (wave blocking in fluid-mechanics 
terminology). The case of waves with negative co- moving frequency has been included 
throughout; these waves are crucial for the Hawking effect and they have been neglected 
in the fluid-mechanics literature on wave blocking. The Hawking effect is a remarkable 
process in which an incident wave generates a wave with negative co-moving frequency, 
with a resulting amplification of the incident wave (this implies an extraction of energy 
from the flow). It has been shown that this process is robust in the presence of 
dispersion (TUl [HI 1121 [131 HH [ISl UHl [IZ] and the linear theory of surface waves falls 
into the class of systems that exhibit the Hawking effect [H]. Experimental evidence 
of the generation of waves with negative co-moving frequency was reported in [22] and 
further experiments are planned. 

We presented analytical results for the deep water/short wavelength case kh ^ 1 
that are more comprehensive than those given elsewhere. A similarity of the state space 
of the waves to that of a thermodynamic system was pointed out; the curves in the state 
space representing the horizon lines are analogous to curves separating thermodynamic 
phases and there is even an analogue of a thermodynamic critical point. 

Schiitzhold and Unruh showed that the regime of gravity-wave propagation in an 
effective Schwarzschild-like metric corresponds to the shallow water limit kh <^ 1 [TS] . 
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The interaction with the white-hole horizon necessarily tunes a wave out of this regime 
{kh <^ 1) into the kh ^ 1 regime in which dispersion causes the effective-metric 
description to break down. In other words, for pure gravity waves the white horizon 
is not dispersive when kh 1 whereas it is dispersive when kh ^ 1. On the other 
hand, waves in the kh ^ 1 regime considered in this paper stay in this regime in the 
interaction with the counter-flow. 

Previous results on the shape of waves at a blocking line (horizon) can be developed 
further for gravity-capillary waves. The horizons for surface waves can be treated as 
examples of saddle-node lines, which also describe caustics in optics. It is well known 
that an Airy function describes both the intensity of light close to an optical caustic of 
the fold type [19], [50] and the water shape for the amphtude of gravity waves close to a 
blocking line [221 EHl E] • It will be shown elsewhere that the Airy function depends on 
a "stopping length" Lg (roughly the width of the arch of the Airy function) which scales 
like Lo ~ qT^/^ (^) . where x* is the position of the white horizon, and f^) , is 

^ \ ax J x=x* ^ ' \ dz / x=x* 

the "surface gravity" at the horizon. An experimental measurement of the Airy shape 
has been carried out by Chawla and Kirby [27J. If surface tension is taken into account, 
we have seen that there is a critical point at the intersection of two saddle-node lines 
(white and blue horizons). We anticipate that the wave at the critical point will be 
described by a Pearcey catastrophe integral [36l Hi] due to the superposition of two 
Airy catastrophe integrals for the two saddle-node lines. 
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